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Abstract 

We study the avoided crossings in the dynamics of quantum con- 
trolled excitations for an interacting two-boson system in an optical lat- 
tice. Specifically, we perform numerical simulations of quantum control in 
this system where driving pulses connect the undriven stationary states 
in a manner characteristic of Stimulated Raman Adiabatic Passage (STI- 
RAP). We demonstrate that the dynamics of such a transition is affected 
by chaos induced avoided crossings, resulting in a loss in coherence of the 
final outcome in the adiabatic limit. 

1 Introduction 

Stimulated Raman Adiabatic Passage, or STIRAP, is a well-known method of 
inducing coherent excitations of quantum systems from the ground state to 
states with higher energy. This is achieved using coherent time-modulated laser 
fields that result in complete population transfer from an initially populated 
ground state to a target state via an intermediate state. STIRAP was first 
proposed by Hioe and coworkers pQ, [2J. A crucial preliminary work by Becker 
et al. [3] achieved efficient vibrational excitation by using a molecular beam as 
an optically pumped active medium and led to the development of the STIRAP 
concept. The theoretical work was formulated by Kuklinsky, Gaubatz, Hioe 
and Bergmann shortly thereafter 0]. STIRAP was further confirmed by ma- 
nipulating the vibrational and rotational degrees of freedom in sodium dimers 
0- 

Since then, STIRAP has been used to control transitions in a diverse range of 
matter-optics systems [6] , ranging from molecular alignment [7] , and molecular 
rotation [8], to the coherent acceleration of ultracold atom systems by trans- 
fer between a stationary and a moving optical lattice [5], and controlled dipole 
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excitations in ultracold bosons subjected to radiation induced double well po- 
tentials [10]. In the latter two cases, STIRAP has also been used to illustrate 
the influence of the underlying classical chaos in the atom dynamics by looking 
at avoided crossings in the Floquet eigenphase spectrum. The exact procedure 
was first described by Na and Reichl [11] for a driven single particle in a box. 

Optical lattice systems have been of great interest in experimental physics 
[12] and [13]) and theoretical physics [14] , [15] for some time. In 1992, Graham, 
Schlautmann and Zoller showed that the center of mass motion of cold atoms 
in an optical lattice could be decoupled from their internal degrees of freedom 
[IB] . Since then, the influence of chaos in optical lattice systems has been an 
important factor in the manipulation of such systems [5] , [17] . Number squeez- 
ing and subpoissonian distribution of atoms in each site in an optical lattice 
have been reported by Itah et al. [18] . The system of interest in our work is an 
optical lattice with a 2-periodic boundary condition that traps two interacting 
bosons. We seek to use this system as a test case for larger many-particle prob- 
lems, where the Hilbert space can be restricted to the subspace of all iV-periodic 
wavefunctions, with large N. Here, we set N to the lowest meaningful value of 
2, since singly periodic wavefunctions do not take into account the probability 
of two particles being in separate wells of the lattice. Furthermore, Fernholtz 
et al. [T5] have recently shown that it is possible to trap a cold atomic system 
on the surface of a toroid and to achieve a two-dimensional periodic potential 
similar to the ringed optical lattice that we shall consider in this work. 

In subsequent sections, we consider the application of STIRAP to an inter- 
acting boson system confined to a one-dimensional optical lattice with periodic 
boundary conditions. A linear version of our system has been implemented 
in experimental studies [20] using ultra-cold atoms. In section [2] we derive 
an expression for the basic model and we discuss the numerical methods used 
to obtain the stationary eigenstates for this system. In section [3] we discuss 
the process by which coherent transitions between selected symmetrized energy 
eigenstates can be achieved for this system. We will also show that avoided 
crossings in the Floquet states can be associated with real transitions of the 
undriven symmetrized eigenstates. Concluding remarks are made in section 3] 



2 The Basic Model 



Our system consists of two atoms (bosons), each of mass to, confined to a 
spatially periodic optical lattice of radius p. The dipole interaction between the 
atom and the optical lattice gives us the atomic Hamiltonian 

where / = mp 2 , Li and 6i are the angular momentum and angle, respectively, 
of the «th particle (i = 1, 2), kq is the lattice amplitude, and uo is the strength 
of the point contact pseudopotential interaction between the bosons. 
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It is useful to write the Hamiltonian in terms of dimensionless parameters 
(L' it #■). For particles interacting with the optical lattice in Eq. (QJ, the angular 
momentum transfer occurs in discrete units AL = 2h. We therefore define 
L\ = fa and fij = 26,. Thus, E{ = K ' = u ' = ^ and f = 4 Wr i, 

where uj r — -Jfj is the recoil frequency. We also scale all other frequencies as 
uj' = . We then drop the primes on the dimensionless parameters and obtain 
the dimensionless Hamiltonian 

H a = L\ + L\ + k [cos6»i + cosfl 2 ] + u 5(6x - 6 2 ). (2) 

We use as a basis set, the eigenstates of the angular momentum operator L\n) = 
n\n) or (#|ri) = "yf^ 6 ™ 9 where integers n range over the values — oo<7i<oo . 

We numerically diagonalize the Hamiltonian in Eq. ([2|) using a nonadap- 
tive finite element method. The 2-particle boson states are obtained by taking 
symmetrized products of single particle states: 

(0i,e 2 \n 1: n 2 ) {s) = -L^imX^K) + (flilnaX&lm)]. (3) 

These states are then used to create a Hamiltonian matrix from Eq. @. The 
eigenvalues E a and eigenvectors \E a ) of the Hamiltonian matrix were deter- 
mined numerically using the appropriate subroutine for diagonalizing real sym- 
metric matrices in the GNU Scientific Library [3T]. Figure [TJ a shows the en- 
ergy levels E\ through £^6- Figure [T]b gives of magnified view of these levels. 
Figures HJa through Oe are the probability distribution plots for the states 
(6 1 ,e 2 \E 1 ), (6 1 ,e 2 \E 3 ), (e^d^), {6 u 6 2 \E 5 ), and (6 1 ,e 2 \E 15 ) respectively at 
the value Ko = 7.287781, the lattice amplitude we use in subsequent sections. 
We chose these states because they have the largest coupling {E a i\ cos(9i)\E a ) 
and lead to robust STIRAP transitions. 



3 Induced Transitions in the Interacting System 

Our goal will be to transition the bosons from the ground state \Ei) of the 
optical lattice to the state I-E15), using \Ei) as the intermediate state. We will 
accomplish this by perturbing the system with two Gaussian shaped radiation 
pulses, the first pulse with carrier frequency Uf = E 15 — E4 and the second 
pulse with carrier frequency w s = E4 — E\. In the presence of these pulses, the 
Hamiltonian takes the form 

H(t)=Ll+Ll + K(t)[cose 1 +cose 2 ]+u S(e 1 -e 2 ), (4) 

where 

n{t) = k + Xf(t) cos(w ft) + X s (t) cos(oj s t). (5) 

Here, the Gaussian amplitudes A/(i) = Aoexp[— (t ~ tj) 2 /4tj\ and X s (t) — 
\oexp[-(t-t°) 2 /4t 2 d ]. 
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As shown in references jS], [5], [TU] and [TT], it is possible to use Floquet 
theory to analyze the effect of the radiation pulses on the system as they pass 
through the system. The only requirement is that the pulse envelopes Xf(t) 
and X s (t) be slowly varying compared to the carrier wave periods 2ir/ujf and 
2tt/uj s , respectively. We can then break the time evolution into N narrow time 
windows, the kth window centered at a time t = th ix . In the kth time window 
we can write the Hamiltonian in the form 



H(t;t k flx ) = L 2 1 +L 2 2 + K(t;t k flx )[cos0 1 +cos0 2 }+uoS(9 1 - 2 ). (6) 



The Hamiltonian H(t; tf ix ) for the feth time window is time-periodic if Uf and 
ui a are commensurate and we can use Floquet theory to analyze the behavior 
of the system in that time window. The value kq — 7.287781 was chosen to 
achieve the commensurability u)f/u> s — 3/2, with OJf = 10.896058668420753. 



The period of the Hamiltonian H (t; t k t ix ) is then given by T = n I jj- + jj- 



The Floquet Hamiltonian, for the fcth time window, is given by Wp{£) = 
Hit; t k fi X ) — i-§i ( m dimensionless units). H F (t) is Hermitian and has eigenvalues 
Q, v and eigenvectors \4> v (t)). The eigenvalues are defined modulo 2tt/T and the 
eigenvectors are time-periodic with period T. The quantum state \ip k (t)}, which 

is a solution to the Schrodinger equation i 9 ^ g ^ = H(t;t k ix )\ip k (t)} can be 
expanded in a spectral decomposition 



where A v gives the contribution of the z^th Floquet state to the evolution of the 
system in a given time- window. The Floquet evolution operator, Uf(T), in the 
basis of symmetrized 2-boson states \n) = \m, n2)^ (see Eq. [3]), is given by 



The Floquet evolution matrix is constructed as follows. To obtain the mth col- 
umn of (n\UF(T)\m), choose initial conditions (m'|^>(0)) = S m ^ m and integrate 
the Schrodinger equation from time t — to time t = T ' . The state at time t = T 
is the mth column of the Floquet evolution matrix. For our system, the integra- 
tion was done using an 8 th order Runge-Kutta Prince Dormand algorithm [22] 
from the GNU Scientific Library [21] . and the Floquet evolution matrix was 
diagonalized using a parallelized LAPACK library through the Scalable Library 
for Eigenvalue Problem Computations (Slepc) [23]- The eigenvalues, e~ invT 
can be used to obtain We obtained the Floquet eigenvalues and eigenvec- 
tors for each time window. The eigenstates in neighboring time windows will 
be approximately orthonormal and we can use this fact to follow each Floquet 
eigenstate and eigenvalue as the pulses move through the system. 

We chose the following parameters for the Gaussian pulses; Ao = 0.2, tf = 
(l/3)i tot , t s — (2/3)t to t, and td = (l/14)i tot . Here t to t defines the total time 
scale for both pulses, and tfi X is expressed in units of ttot unless otherwise stated. 





(7) 



(n|f/ i ,(T)|m)=^e-^ T (n|^(0))(^(0)|m). 



(8) 
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Figure [3] shows the Floquet eigenvalues of the relevant Floquet eigenstates as 
the system evolves in adiabatic time t fi X . The relevant eigenstates are the ones 
isomorphic to the states connected by the STIRAP pulses viz. \Ei), \E±), and 
\Ei5). We notice that the Floquet eigenvalues are degenerate at tfi X =0 and 
tfix = Uot as expected. The Floquet states and corresponding eigenvalues are 
labeled alphabetically as follows: 

1. The cigenphase whose corresponding Floquet eigenstate is supported by 
the undriven state [I-E4) — l-Eis)] at tfi x — is labeled as £Ib and the 
Floquet eigenstate as |^>s). 

2. The eigenphase whose corresponding Floquet eigenstate is supported by 
the undriven state [I-E4) + I-E15)] at tfi X — is labeled as tie and the 
Floquet eigenstate as \4>c)- 

3. The eigenphase whose corresponding Floquet eigenstate is supported by 
the undriven ground state \E±) at tfi x = is labeled as flo and the 
Floquet eigenstate as \4>d)- 

The evolution of the Floquet eigenvalues f2s, f^cs and £Id are shown in Fig. 
[3] as the pulses pass through the system. If the system evolves adiabatically, it 
stays in the state |</>d), but the amount of support it receives from each of the 
states \Ej) will change at each avoided crossing that \4>d) encounters. In Fig. [31 
there appears to be a traditional 3-level avoided crossing at about tti x ~ 0.5 ttot- 
However, there are additional avoided crossings involving that can change 
\4>d)- In Fig. [31 we show the dependence of the three Floquet states \(/>b), |0c) 3 
and \4>d) on the energy eigenstates \Ej) as the pulses pass through the system. 
We see that \4>d) starts out fully supported by \E\) then partially changes its 
support from \E\) to \E\$) and finally, after the pulses have passed, ends with 
equal support from \E\) and \E±). State \4>b), on the other hand, begins with 
equal support from \E\) and \E±) and ends totally supported by I-E15). A truly 
adiabatic evolution of the pulses, which would keep the system in Floquet state 
\4>d) the whole time, would not accomplish our goal of transitioning the system 
from \E\) to I-E15). The problem arises because additional avoided crossings 
occur that pull the system off the traditional STIRAP path. The influence of 
an additional small avoided crossing between VLb and flu at about tfi x = 0.55 
is clearly seen. This additional small avoided crossing is a manifestation of 
classical chaos in the quantum dynamics similar to the ones seen for the double 
well system in [TU] . 

The dynamics can be analyzed in more detail by using the Landau-Zener 
formula to calculate the probability of a transition at an avoided crossing. The 
probability P v y for an avoided crossing between two Floquet eigenphases fi„ 
and fi„/ to be crossed is given by [24] , [25] 



2 1 



P u y = exp - 



(9) 
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where 8£l u y is the (minimum) spacing between tt u and Q! u at the avoided 
crossing and Y v y is the magnitude of the diabatic rate of change (slope) of the 
Floquet eigenphases. Thus, 



dt dt 



(10) 



where ^ L is the slope of the eigenphase curve fi„ in the neighborhood of the 
avoided crossing. Equation ([9]) can be simplified to P v y — exp [— ttot'yv,v']i 

where j v y = 2r > t — . 

In order for a crossing to be traversed adiabatically, P UM i « 0, and the 
actual time scale of the STIRAP must be adjusted accordingly. Thus, the 
transfer probability P v y will be very small if t to t > 1 /ivy ■ For the BD avoided 
crossing, we estimate the gap SttsD to be 1.984 x 10~ 3 , and Tbd to be 0.135, 
concluding that t to t > 2.184x 10 4 . For 85 Rb atoms confined in a one-dimensional 
optical lattice by detuning away from the Di transition line, the recoil frequency, 
Lu r is about 24 KHz [20]. The characteristic time scale here is l/(4cj r ), or 
1.03 x 10 -5 seconds. We can plug this value to the minimum value(s) of t to t to 
get the actual time. Thus, for the 85 Rb atom, we get t to t > 0.225 sec for the BD 
crossing. If the time scale for the stirap is faster than 0.225 seconds, then the 
crossing is traversed diabatically and traditional 3-level STIRAP will be seen. 
For 87 Rb atoms in the F — 2, rap = 2 state confined in the toroidal magnetic 
trap of Fernholtz et al, the resonance condition, as well as one dimensional 
confinement, are met at a radius ~ fxm [19] . and can be adjusted to replicate 
the time scales of the optical lattice. 

We now compare the results obtained from Floquet theory to numerical 
simulations of the actual dynamics of the system as it evolves through time. The 
full multilevel Schrodinger equation for this system, starting from the undriven 
ground state IV'(O)) = \E\) , can be solved numerically for the system as it evolves 
from t = to t = t to t ■ The rate at which the avoided crossings are traversed can 
be controlled by controlling t to t- Figures [5]a through [5Jd show the numerical 
time evolution of the wavefunction of the two-boson system \ip(t)}, starting from 
the ground state For small values of ttot: the system traverses the BD 

avoided crossing diabatically, and a near-complete population transfer to I-E15) 
occurs due to the preceding 3-level avoided crossing, replicating the traditional 
STIRAP process. This can be seen in Figure [SJa, where t to t = 7200. As we 
increase ttot towards (and beyond) 2.184 x 10 4 , the time evolution approaches 
that of the Floquet state \4>d) as shown in Fig[4j Figure[5]d shows \ip(t)) as it 
evolves in time for t to t = 90, 000. The centroids of the components (Ei\tp(t)) are 
identical to the components (i?i|</)£>) in Fig0J The oscillations in the probability 
occur due to nonadiabatic effects as demonstrated by Berry [25], and decrease 
in amplitude as we move further into the adiabatic regime. 
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4 Conclusion 



We have analyzed the dynamics of interacting two-boson systems for ultracold 
alkali metal atoms in an optical lattice with periodic boundary conditions. We 
have demonstrated the feasibility of a controlled excitation of the system into a 
higher energy state using STIRAP, induced by time dependent modulations of 
the optical lattice. 

For sufficiently large amplitude modulations, the effects of the underlying 
classical dynamics [27] start to manifest themselves through small avoided cross- 
ings between the involved Floquet eigenphases. The STIRAP pulses were tuned 
to connect very high energy states (the final state being the fifteenth energy 
level). Avoided crossings between the other Floquet states connected cause the 
outcome to differ from traditional three level STIRAP. By traversing these ad- 
ditional small avoided crossings diabatically, in order to eliminate their effect on 
the system, we obtain the outcome expected for a three-level STIRAP process. 
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1 Energy curves of the first 66 lowest-energy states of both even 
and odd parities for the two-boson system, (a) The energy levels 
plotted as a function of kq for interaction amplitude uq = 23.0. 
The region of interest is highlighted by a box. (b) Magnified view 
of the region of interest boxed in Fig. l.a. The levels being con- 
nected by STIRAP for this particular value of kq are indicated. 
The value of ko has been adjusted so that — £ = | 

2 Wavefunction plots for states (a) |£x), (h)°\E 3 ), (c) \E 4 ), (d) 
\E 5 ), (e) \E 15 ) 

3 Floquet quasienergies f2„ as a function of adiabatic time tfi X /t to t 
for A = 0.2. The labels Qb, ftc, and J1d denote the floquet 
quasienergy curves corresponding to the Floquet eigenstates 10b), 
1 0c), and \4>d), respectively 

4 Plot of the Floquet eigenfunctions in the undriven Hamiltonian 
representation ((E^b-d))- The components of the Floquet 
states in each energy level are numbered. Note the influence 
of the avoided crossings in \4>d) 

5 Plots of |(i£j-|^(t))|| 2 as a function of time of the system as it 
evolves under the influence of STIRAP pulses from the ground 
state IV'CO))! = | -Ei) for different values of ttot- All units are 
dimensionless. (a) t to t = 7200. (b) t to t = 12,000. (c) t to t = 
24,000. (d) t tot = 90,000 
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Figure 1: Energy curves of the first 66 lowest-energy states of both even and odd 
parities for the two-boson system, (a) The energy levels plotted as a function 
of kq for interaction amplitude uq — 23.0. The region of interest is highlighted 
by a box. (b) Magnified view of the region of interest boxed in Fig. l.a. The 
levels being connected by STIRAP for this particular value of k are indicated. 
The value of kq has been adjusted so that ^ = |. 
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Figure 2: Wavefunction plots for states (a) \Ex), (b) \E 3 ), (c) |£ 4 ), (d) \E 5 ), (e) 
\Ei 5 ). 
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Figure 3: Floquet quasienergies VL V as a function of adiabatic time tfi X /t to t for 
A = 0.2. The labels fls, and £Id denote the floquet quasienergy curves 
corresponding to the Floquet eigenstates \(/>b), \<t>c), an d \4>d), respectively. 
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Figure 4: Plot of the Floquet eigenfunctions in the undriven Hamiltonian rep- 
resentation ((Ei\(j)B-D))- The components of the Floquet states in each energy 
level are numbered. Note the influence of the avoided crossings in \<Pd)- 
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Figure 5: Plots of |(.Ej|^(i))|| 2 as a function of time of the system as it evolves 
under the influence of STIRAP pulses from the ground state IV^O))! = \E\) 
for different values of t to t- All units are dimensionless. (a) t to t = 7200. (b) 
t tot = 12, 000. (c) t tot = 24, 000. (d) t tot = 90, 000. 
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